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I. INTRODUCTION 


The detection of the Higgs particles is one of the most important objectives of the Large 
Hadron Collider (LHC). Charged Higgs bosons are predicted in extended versions of the 
Standard model (SM), like two-Higgs-doublet models (2HDM) and the Minimal Supersym¬ 
metric Standard Model (MSSM). Discovery of such an additional charged Higgs boson will 
immediately indicate physics beyond the SM, unlike the case of the neutral Higgs boson. 
Hence, there is strong theoretical and experimental motivation for exploring the mechanisms 
of the charged Higgs boson production. 

The charged Higgs boson H could appear as the decay product of primarily produced 
top quarks if the mass of H ± is smaller than rn t — m&. For heavier H ± , the direct H ± 
production mechanisms at hadron colliders have been extensively investigated. At the LHC, 
the primary charged Higgs boson production channel is gb —> H~t [1] *. The study [3,4] 
shows that this production mechanism can be used to explore the parameter space of MSSM 
for rn H ± up to 1 TeV and tan/? down to at least ~ 3, and potentially to ~ 1.5. Therefore, 
it is necessary to calculate and implement also the loop contributions to gb —» H~t for more 
accurate theoretical predictions. 

The production mechanism of the heavy charged Higgs boson with heavy top quark 
has been studied two decades ago [5]. In order to resum the possible large terms like 
log (Q 2 /ml), the bottom quark parton distribution function (PDF) is introduced. Although, 
there are some doubts about the bottom parton description [6], the detailed study on the 
Pr.b distribution [7] argues that the bottom parton description is reliable for the inclusive 
tH~ production process at LHC. 

In earlier literature, the contribution of the initial-gluon process gg —> H~tb, which is 
only part of the next-to-leading order (NLO) QCD corrections to gb —> H~t, has been 
calculated [8]. The supersymmetric electroweak corrections arising from the quantum ef¬ 
fects which are induced by potentially large Yukawa couplings from the Higgs sector and 
the chargino-top(bottom)-sbottom(stop) couplings, neutralino- top(bottom)-stop(sbottom) 


*See Ref. [2] for the discussion on other charged Higgs boson production mechanisms. 
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couplings and charged Higgs-stop-sbottom couplings are also studied [1,9], which can give 
rise to a 15% reduction of the lowest-order result. In Ref. [10], the electro-weak corrections 
to the process are also discussed. After the submission of this paper, another calculation 
on the QCD correction appeared [7], in which the corresponding results seemed compatible 
with ours. The study on the SUSY-QCD effect for this process is also done [7,11]. In this 
paper, we would present the detailed and complete NLO QCD corrections to gb —> H~t. 
We should note here that the results presented in this paper are for the process bg —> tH~\ 
they are the same for the charge conjugate process bg —> H + t. 

The arrangement of this paper is as follows. Section II contains the analytic results, and 
in Section III we present numerical examples and discuss the implications of our results. 
The lengthy expressions of the form factors are collected in the Appendix. 


II. ANALYTIC EXPRESSIONS 

Including the NLO QCD corrections, the cross sections for PP —» tH~X at the CERN 
LHC can be written as 

a = a LO + a Vir + a Real , ( 1 ) 

where o LO is the cross section at leading order (LO), <J Vir and a Real are cross sections from 
NLO QCD corrections arising from virtual and real processes. 

A. LO cross section 

The Feynman diagrams for the charged Higgs boson production via b(pi)g(p 2 ) —> 
t(ki)H~(k 2 ) at the LO are shown in Fig.l. The amplitudes are created by use of Feynarts 
[13] and are handled with the help of FeynCalc [14]. As usual, we define the Mandelstam 
variables as 


s — (p\ + P 2 Y — {k\ + k 2 ) 2 , 
t = (pi - h) 2 = (p 2 - k 2 ) 2 , 

11 - (p\ - h) 2 = (p> - ki ) 2 . (2) 
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The amplitude at the LO could be written as 

6 

MlO — E ^ \ c l^2i-l + C 2 M 2 j], 

i =1 

where the non-vanishing form factors are 

1 1 

^2 o 5 

mf — u s 

2 

^3 2 7 

— u 

2 

^5 =-, 

S 

where Ci = gmg ^% mttan/9 and c 2 = gmg ^ mt cot ‘ . In Eq. (3) M 2i and M 2l+X (i 
standard matrix elements which are defined as 

Mi = u(ki) /:(p 2 )Pru(p 1 ), 

M 2 = u(ki) 4(p2)P L u(pi), 

M 3 = m(/ci) i> 2 fc{k 2 )P R u(pi), 

M 4 = u(Aii) i> 2 /(/c 2 )P l m(p 1 ), 

M 5 = • e(p 2 )| 

M 6 = it(fci)Pi,tt(pi)fci ■ e{p 2 ), 

My = m(/Ci) £) 2 P R u(p\)k\ ■ s(p 2 ), 

M s = u(k i) jt> 2 P L u(pi)ki ■ e(p 2 ), 

Mg = 'u(/h)P R 'u(pi)pi • e(p 2 ), 

Mio = m(A;i)Plm(pi)pi • £(^ 2 ), 

Mn = h(Pi) £) 2 Pru{pi)pi • e(p 2 ), 

M i2 = w(/ci) P 2 Plu(pi)pi ■ s(p2), (5) 

where the color matrix T a has been suppressed. In this paper, we perform the calculations in 
Feynman gauge and in d time-space dimensions with d — 4 — 2e. For simplicity, throughout 
the paper we omit the bottom quark dynamical mass but keep the mass-term only in Yukawa 
couplings. The discussion on the error induced by massless bottom quark approximation 
can be found in Ref. [12]. 


( 4 ) 

= 1-6) are the 
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The LO cross section can then be written as 


with 


d(JLO 

dx\dx2 


dcr 0 Gb/A(xi, fif)G g /B(x2, df) + [A 


B\, 


( 6 ) 


dcr° 


111 
24 4(1 — e) 2s 


\M LO \ 2 d<& 2 


( 7 ) 


where the factor and are the color and spin average, respectively, and two-body 

phase space is 


d <&2 


1 /4 7T\ e 1 

87t V s / T(1 — e) 


A 




1/2—e 

n“ e (l 


n) e dn, 


( 8 ) 


with v = |(l + cos$), where 0 is the center-of-mass scattering angle between p\ and k\. Here 
A is the two-body phase space function 


A (a;, y, z) = x 2 + y 2 + z 2 — 2 xy — 2 xz — 2 yz. 


(9) 


B. Virtual corrections 


The Feynman diagrams for the NLO virtual corrections are shown in Fig. 2. The virtual 
diagrams with the self-energy insertions on the external legs are not shown. They can be 
obtained by inserting diagram (a)-(d) of Fig. 2 into the external legs of the LO diagrams 
in Fig. 1. At the same time, the diagrams containing counter-terms can be easily got by 
inserting corresponding counter-terms into external legs, internal propagators and vertex of 
the LO diagrams. 

In order to remove the UV divergences, we have to renormalize the strong coupling 
constant, the Yukawa coupling constants, the quark masses and the wave functions of quarks 
and gluon. The strong coupling constant and the gluon wave function are renormalized in 
the MS scheme as 


5g s 


9s 

z 9 


cv 

--^(2C a -(3 0 )A, 


( 10 ) 
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with A = i — 7 ^ + log( 47 r), /3 0 = (11 C A — 2n/)/3, CA = 3 and Cp = 4/3. Here, rif = 6 is 
the number of fermions. 


We renormalize wave function and mass of the bottom quark in the MS scheme 


Sniij 

rn h 


Z b 


a s 

47T 


3C F A, 


->*■ 


in) 


At the same time, we will renormalize the top quark mass in two schemes: on-mass-shell 
(OS) and MS. 


Snit 

m t 

Sm t 

m t 


-r? c r 


A + - - log(m t 2 //i 2 ) 


a, 


- 3CpA, in MS scheme. 

47T 


in OS scheme 


( 12 ) 


Hereafter we will refer OS and MS schemes to the different mass renormalization of the top 
quark. The comparison between the results in these two schemes will be discussed in the 
numerical section. The wave function renormalization constant of top quark is chosen the 
same with that of b quark 


£t = 


a s 

47T 


C F A. 


(13) 


The renormalized virtual amplitude can then be written in the following way, 


M ren = M a + M b , (14) 

where M A is the amplitude from the diagrams (e)-(o) of Fig. 2 and Mb is the amplitude 
from the diagrams which contain self-energy insertion on the external legs or counter-terms. 
The Mj 4 can be written as 

Ma = ±M\, (15) 

i=e 

where i represents the diagram index of Fig. 2. For each diagram i, we can generally write 
the amplitude as 

6 

M\ = ]T fAdM^ + c 2 M 2l ], (16) 

3 =1 
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where the non-vanishing form factors fj are given explicitly in Appendix and the Mj are the 
standard matrix elements given in the previous subsection. 

The Mb can be written as 

M b = M l B + M %, (17) 

where 



Ml = + c 2 f?M 2i \. 


Here the non-vanishing form factors — l — 12) are 

,1 _ f2 _ Sm t 

J S J S 9 5 


mf — u 

—2 mtSmt 8m b 


(rrif — u) 2 


rrib I rrif — u 
5rn t T 1 


4 _ -2 m t 5m t 5m t 

J S / 9 \ 9 ' 9 

m t \_mf 


fs= 2 


/s=2 


—2 m t 8m t 5m b 1 

(m 2 — w) 2 nib m t ~ u 

—2 m t Sm t Sm t 1 

(m 2 — w) 2 m t m 2 — u 


ft = — 


A 10 = - 


s m b 
2 8m t 


After squaring the renormalized amplitude and performing the spin and color summa¬ 
tions, the partonic cross section with virtual corrections can be written as 

da Vir _ J_ , , r . _ 


dxidx 2 


— 2 Re ^2 (M+ n M LO ) l J -f) + [A <->• H], 


where A, H denote the incoming hadrons. 

After the renormalization procedure described above, da Vir is UV-hnite. Nevertheless, 
it contains still the soft and collincar divergences. The soft divergences will cancel against 
the contributions from real-gluon radiation [see Eq. (24)]. The remaining collinear diver¬ 
gences in the real-gluon-emission processes will be removed by the redefinition of the parton 
distribution functions (PDF) (mass factorization) [see Eq. (27)]. 
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C. Real corrections 


There are three kinds of real corrections to the processes bg —> tH~: gluon-radiation 
[bg —> tH~g], initial-gluon [ gg —> tH~b ] and initial-active-quark [bq(q) —> tH~q(q) and 
gg —> tH~b, where q stand for the active quarks which are treated as light for PDF evolution, 
in practice the light quarks other than u, d, s can be omitted due to the low luminosity]. All 
real corrections are related to the 2 —> 3 processes. 

In this paper, the 2 —> 3 processes have been treated using the two cut-off phase space 
slicing method (TCPSSM) [16]. The method is briefly described in the following. Two 
small artificial constants 8 S , 8 C are introduced, and the three-body phase space can firstly be 
divided into soft and hard regions according to whether the gluon energy is less than 8 s \/s/2. 
Secondly the hard region is further divided into hard collinear and hard non-collinear regions 
according to whether the magnitude of Pi-Pj is less than 8 c s/ 2 (jg , p 3 are the possible collinear 
momenta). In the soft and collinear regions, the phase space integration can be performed 
analytically in d — 4 — 2e dimensions. At the same time, in the hard non-collinear region, 
the phase space integration can be calculated in four dimensions by standard Monte Carlo 
packages because the integration contains no divergences. Obviously, the final physical 
results should be independent on these artificial parameters 8 S and 8 C1 which offers a crucial 
way to check our results. Therefore, the real corrections can be written technically as, 
according to the phase space slicing, 


o R eai = o s + a CoU + a fin 


( 22 ) 


with 


a Coll = a Coll + a fac } (23) 

where a s , a Co11 and <r-' m are cross sections for the direct calculations of the 2 —» 3 processes 
in soft, hard collinear and hard non-collinear regions, and the a^ ac is the counter-term from 
factorization procedure. After adding virtual contributions and a s , the double poles are 
canceled, and the remaining singularities together with a Co11 are canceled by a* ac . In the 
following, a s , a Co11 and (7^ n will be presented respectively. 
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1. Cross section in soft region 


The Feynman diagrams of the gluon-radiation process bg —> tH~g are shown in Fig. 3. 
Diagrams except (c) and (d) contribute to the cross section in the soft region. We may write 
the cross section as 


dcr s 

dxidx 2 


— d^Gb/A^xi, g,f)Gg/s{x2i /if) + [A <-> B ), 


a s = a 


a s T(1 — e) / 4:H(i 
2^T(1 -2e) I" 


2 \ e 


A 1 as 

7X + — + Ao1 ’ 


(24) 


where 


^2 = -- 


rnf — t 


+ 12 - 


rnf — u 


3 (Ei — (3 cos 9)s (Ei + /3 cos 9)s 
4 m 


+ 12 , 


16 4mf 
A{ = 1 


+ 


t 


3 (El -p 2 )s 3(E V -p cos 9)s 

,„ 2 

+ 12 


(Cl -:2lOg ——) 


mf — u , ^ , 2 . , 2 

(C 2 - 2 log —-=) + 24 log 


(E\ + /3 cos 9)s 
1 , 2 4 4 16 4 ml 

A °~2 l0g ^ 2 “ l0g ^ + T (E 2 - p 2 )s 


- log ~A s 2 , 
s s 


mf — t 


3 (Ei — (3 cos 9)s 


-12 
Ci = log 
C 2 = log 


mf — u 


(Ei + (3 cos 9)s 
(Ei — (3 cos 9) 2 
E 2 - (3 2 ’ 

(Ei + (3 cos 9) 2 
E 2 - (3 2 ’ 

.2 El-(3 


(C 3 — 2 log 2 
(C 4 — 2 log 2 


mi ^ Bp Cx + /5 
(2 ‘° g + 8 
2 . 




+ 2Ci log 


2 


+ 2C 2 log -—-=) + 24 log 


/-/ i 2 I 1 ] 2 El + [3 07 ' / 

C3 = -‘° g g.-ffcos+ F 08 £T+" 2 ' !2( 

C1 ~ ” ‘° g £ 1+ ^cos« + 2 l0g “ 2( 

P = *3(1 ,mt/s,m 2 H ±/s), 


/3cos9 + /3: , /?cos0 + /3 

) + 2 li 2 ( 


Ei-(3 

p COS 0 + P 


Ei — P 


) + 2 li 2 (- 


Ei — P cos 0 7 ’ 
~P cos 6 + /+ 

Ei + P cos 0 


Ei = \P 2 + 


2 m? 


( 25 ) 


with cos 9 defined in 11-A. 



2. Cross section in hard collinear region 


As discussed above, the real corrections in collinear region <j Co11 contain divergences 
because the partons are massless. In order to remove these kinds of collinear singularities, 
we introduce scale dependent parton distribution functions in MS convention 


G c /p(x, ixf) — G c /p(x ) — 


CK S r(l - e) 

/4vr/i 2 yi 

_27t T(1 — 2e) 

V y} ). 


r4 dz 


lx Z 


Pc C f(z)G c >/p(x/z), 


(26) 


where P cc '(z) are splitting functions and c[c r ] represent b or g [partons which can split into 
c]. After substituting the G c /p{x) in the lowest order expressions by G c /p(x, /if), we can 


obtain <jf ac . Adding a Co11 and a* ac , we can get the final results for a Co11 
a s T(1 -e) (4 :-k n 2 \ " 


da Co11 

dx\dx2 


27T T(1 — 2e) y s 

w , A T( b -»■ b 9 ) , A T(g -»• 99 ) 

X | I 


\^Gb/A{Xl, i7f)Gg/B(x 2 , /i/) 

+ A s 0 c (b - bg) + A£{g - gg) 


+ Gb/A^Xi , Hf)Gg/B(x 2 , Hf) + (j&/a(^1, t L f)Gg/B(x2, /i/)| do 

+ [A ^ 5], 


where [16] 


Af(6 


• 4 “ = - 4 “ ln (^J 

bg) = Cf( 2 In + 3/2) 


gg) = 2N In 6 S + (HAf — 2rif)/6 

with rif — 5 is the number of massless fermions. Here 

G c/ B,A( X i^f) I CC — G c'/B,A( X /y^f)Pcc'(y) 

r> Jx y 


with 


PM = p j(y) ln [£ —--t) _ p ij(y )* 
v y pf j 


where 


P„(z) = C, 


1 + z 2 

' 1 - Z 




(27) 


(28) 

(29) 

(30) 


(31) 


(32) 


9 


(33) 

(34) 



with N = 3. 





C, 


1 + (1 -zf 


- C F z 
2N 


z 1 — z 

+- +z{ 1 


1- z z 


0 

\ P + (1 - *) 2 ' 

-z(l - z ) , 



(35) 

(36) 

(37) 

(38) 

(39) 

(40) 


3. Cross section in hard non-collinear region 

As described above, the cross section in hard non-collinear region can be easily 
obtained by Monte Carlo phase space integration in four dimension. It can be written as 

drrf in 

-—-— = X/ G c /a( x i, d'f)G c '/ B ( x 2, ^f)\cc' —> tH X | 2 d<h 3 + [A <-»• B ], (41) 

ClX\CtX 2 ( , (J 

where c, d run through gluon and light quarks and the three-body phase space < 1>3 is within 
the hard non-collinear region. In this paper, all Monte Carlo phase space integrations are 
performed by package BASES [17]. 


III. NUMERICAL RESULTS AND DISCUSSION 

Our numerical results are obtained using CTEQ5M (CTEQ5L) PDF [18] and 2-loop (1- 
loop) evolution of a s (fi ) for NLO (LO) cross section calculations with A^ = 226 (146) MeV. 
As we have mentioned in the second section, the dynamical mass of the bottom quark has 
been set to zero, which means that we treat the bottom quark as usual massless partons. 
In the MS scheme, 2-loop evolution of the quark masses is adopted, and the pole masses 
of bottom and top quarks are taken as 4.7 GeV and 175 GeV. In the OS scheme (defined 
in this paper), the top quark mass is equal to the pole mass and the bottom quark mass is 
the same with that in the MS scheme. For simplicity, the renormalization and factorization 
scales are taken to be the same. 


10 



In Fig. 4, we show the cross sections as a function of 5 S with 5 C = <5 s /50 in the MS scheme. 
It is obviously true that, as stated in II.C, the physical cross sections are independent on 
the artificial parameters 8 S and 5 C . 


A. Theoretical uncertainties 


There are many theoretical uncertainties in the calculation of the cross sections, for 
examples unknown scales of renormalization and factorization as well as the different choice 
of renormalization schemes. We can define two quantities 5 and A to measure such kinds of 
uncertainties as 


W = 

A(/i) 


°~Os(/^) °Msi[ u ) 

<7 os{f 1 ) + a Ws(^) 
= u(2p.) - o-(/V2) 
cr(2/i) + a(n/2)' 


(42) 

(43) 


It is known that in perturbative calculations, the physical results are independent on 


the renormaliztion schemes provided that we can expand the perturbative series to infinity. 
Otherwise, the results do depend on the renormalization schemes. Therefore 6 is due to 


the unknown parts of higher order effect of the perturbative expansion series, and it can be 
enhanced by the large term log(/x 2 /m 2 ), especially for large tan/3 case, where the yukawa 
coupling containing rn b dominate the contribution. For such case, <5 can be the measure for 
the entire theoretical uncertainties. Explicitly, for large tan f3 limit 


°lo « rn 2 b , 

°lo oc m 6 (/i) 2 . (44) 


Therefore 


ml - m fe (/i ) 2 ^ A 
ml + m b (ii) 2 2 


(45) 


with (at one-loop order) 


A 


2—(log 

7T 




(46) 


which is the quantity entering also the relation between the MS quark mass and the corre¬ 


sponding pole mass, 
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(47) 


m(/i) 2 = m 2 [1 — A]. 

For /j, = fi 0 = m,H± + m t , 5 lo — 30% if we choose m#± = 200 GeV. 

At NLO, we can write 

a NLO oc m 2 (i — A + 5), 

°NLO oc m b (n) 2 (l + B), (48) 

where B is the 0(a s ) radiative correction to the LO cross section in the MS scheme. Hence, 
one finds 

A B 

Snlo — • (49) 

For the same ju in the last paragraph, we have B ~ 0.3, and 5nlo — 10%. 

For the small tan /3 case, where the yukawa coupling containing top quark mass dominates 
the contribution, 5 can be obtained in the similar method as above. For n — n$ — 375 GeV, 
5lo — 10% and Snlo — 3%. For the intermediate value of tan/3, 5 lies between the large 
and small tan /3 case. 

From above discussion on S, we can see that 5 at NLO is not large for the small tan f3 case. 
Therefore the theoretical uncertainties from varying the renormalization and factorization 
scales might be more important. In Fig. 5, we show in (a): the LO and NLO cross sections 
in OS and MS schemes and in (b): the relative deviation S as a function of renormalization 
and factorization scales fi/no for m#± = 200 GeV and tan/3 = 2. From the figure we can 
see that the direct calculation confirms the above estimation for <5. Furthermore, the NLO 
results reduce the scale-dependence in both schemes. From the figure we can also calculate 
Alo,nlo — 6%, 0.2% in OS scheme and A LO:NLO — 17%, 5% in MS scheme for /j = /j 0 GeV. 

From the discussion on 5 and A, which act as the measures of the theoretical uncertain¬ 
ties, we can see that the theoretical uncertainties of the cross sections at NLO are much 
smaller than that at LO. 


B. Numerical results in the MS scheme 

In this section, we will give the numerical results in the MS scheme, in which the 
large terms like a s log(m 2 //r 2 ) have been resummed into the running of the b quark mass. 
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Therefore it gives usually more stable results with respect to missing higher-order terms. 


In Fig. 6 we show the K-factor, which is defined as 

&nlo 


K = 


&LO 


(50) 


as a function of the charged Higgs mass with the renormalization and factorization scales 


H — no- In the MS scheme, if the identical Yukawa couplings are used at LO and NLO 
calculations, the K-factor does not depend on tan/?, which is not true in the OS scheme 
because the different mass renormalization constants of top and bottom quarks spoil it. In 
Fig. 6, the different contributions to K-factor from improved Born (which is equal to 1 if the 
difference between the LO and NLO PDF and a s (n ) is omitted), virtnal+gluon-radiation, 
initial-gluon, bq (q) (q stand for the light quarks) and qq are also shown. From the figure, 
we can see that K-factor from improved Born contribution is around 1.2 ~ 1.3. K-factor 
from virtual+gluon-radiation contribution is from 0.7 to 0.9 when the charged Higgs boson 
mass varies from 200 GeV to 1000 GeV. The initial-gluon and bq(q) contributions to the 
K-factor are negative, and they vary from ~ —27% to ~ —24% and ~ —5% to ~ —14% 
respectively. The qq contribution to the K-factor can be neglected, the magnitude of which 
is smaller than 3% for all charged Higgs boson mass. Adding all the contributions, we can 
see that the K-factor varies from ~ 1.6 to ~ 1.8 when charged Higgs mass increases from 
200 GeV to 1000 GeV. 

In Eq. (50), if we replace the two-loop evolution quark masses by one-loop evolution one 


in the calculation of the LO cross section, in the MS scheme, the K' can be expressed as 

+ m b,(2 ) tan4 P 


K' = K- 


(51) 


m t,(l) + m b,( 1) tan4 P 

where the subscripts of quark masses are which kind of quark mass evolutions are used. In 
Fig. 7, the K' is shown as a function of charged Higgs boson mass with tan/3 = 2, 5,10, 30, 
respectively. The dependence on tan/3 can be explained by Eq. (51). 


To summarize, the next-to-leading order QCD corrections to charged Higgs boson asso¬ 
ciated production with top quark through bg —> tH~ at the CERN Large Hadron Collider 
are calculated in the minimal supersymmetric standard model and two-Higgs-doublet model 
in the MS scheme. It should be noted that in MSSM, the SLISY-QCD corrections arising 
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from the virtual gluino and squarks should be also included. For some specific parameters, 
the SUSY-QCD can be significant [7,11]. From the calculations, we can see that the NLO 
QCD corrections can reduce the scale dependence of the LO cross section. The K-factor 
does not depend on tan (3 if the same quark running masses are used in the NLO and LO 
cross sections, and varies from ~ 1.6 to ~ 1.8 when charged Higgs mass increases from 200 
GeV to 1000 GeV. 
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V. APPENDIX 


In this appendix, we will give the non-vanishing form-factors in Eq. (16). For complete¬ 
ness, we give firstly the definition of loop integrals and its Lorentz decomposition: 

(2vr/i) 4 - d r 1 




in- 


d a q 


QoQi 


(52) 


C 0 ^-MplPi2,pl,ml,rnl,rnl) = 


(2vr/i) 


4 -d 


m- 


>d Qi-iQv 


(53) 


Do-^u{ pI, Pi2, P 32 ,pl,pl,P 3 i, rnl, ml, m 2 2 , m 2 3 ) 


( 2?r p) 4 d f d d !; 

m 2 J q Q 0 QiQ2Q3 


(54) 


with 


Qo = q 2 - ml + ie, Qi = (q + Pi ) 2 - m 2 + ie, p i;j = ( Pi - pj ) 2 (55) 

and 
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C/J, — Pl/xCi + P2[iC 2 

C IW = g^Coo + Pl^PluCn + ( Pl^u + P2,xPlu)C 12 + P2»P2vC 22 

Dfj, = Pi^Di + p 2 fl D 2 + p 3 l _ L D 3 

D IIU = g, lv D m + p ltl p lu D u + {pi^p 2u + P2hP\v)D\2 + {pi/Jhu + p 3 »Piv)D 13 

J rp 2 ^p 2 v D 22 + (p 2 ll p 3l , + p 3 f iP 2 u)D 23 + p 3 / _iP 3 u D 33 . (56) 

For simplicity, we define abbreviation for B l Q {i = 1 — 7), C l x (i = 1 — 8) [x stands for the 
subscript defined in Eq. (56)], D l 0 (i = 1 — 3) as 

B\ = B 0 (0, 0, m t 2 ), 

Bo = B 0 (0,m 2 t ,m 2 t ), 

Bo = B 0 (m 2 H ±,0,mt), 

Bo = Bo(rnf , 0 , ml), 

Bo = B 0 (s, 0, 0), 

Bo = B o (t,0, rrif ), 

B 7 0 = B 0 (u,0,m 2 t ), 

Cl = C x ( 0,0, 5, 0,0,0), 

Cl = c x (0, m 2 H ±, t, m 2 t ,m 2 t , 0), 

C't = C x (m 2 H ±, 0, t, m 2 t , 0, 0), 

C A X = C x (m 2 H ±, 0, u, ml, 0, 0), 

C b x = C x (m 2 H ±,ml, s, 0, m 2 t , 0), 

Cl = C x (ml, 0, t, ml, 0, 0), 

C 7 X = CJrnl, 0, u, 0, m 2 , m 2 t ), 

Cl = C x (ml, 0, u, ml, 0, 0), 

Dl = D 0 (m 2 H± , 0, ml, 0, t, u, 0, m 2 t , m 2 t , 0), 

D'l = Do(m 2 H ±,m 2 t ,0,0,s,t,0,m 2 t ,0,0), 

Dq = D 0 (m 2 H ±, ml, 0, 0, s, u, 0, ml, 0, 0). (57) 

For diagram (e) in Fig. 2, we can write the form factor as 
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(58) 



91 = m«f>(C 0 5 (l - e) + Cf - e(C« + C|)), 

92 = | = Co( m //± - 2m?) - B 0 * - B§ - m?(C? + C?) + e(Bj + S C|). 

For diagram (h) in Fig. 2, we can write the form factor as 

f. = Cf—A - -g, 

87 T z {mf — u ) 

9i — m t (—1 + e )( m t ~~ u )(Cq + C*i + C 2 )* 

£2 = y = -C'ot m l± + (- 2 + e )' m t ~ eu ] + 5 o 
+ m t (C1 + ^2) — e [-^o + ( m t — ' u )(C < i + ^2)]- 


(59) 

(60) 

(61) 

(62) 

(63) 

(64) 

(65) 


For diagram (i) in Fig. 2, we can write the form factor as 


167r 2 (m| — u)' 


gi = 3 m t {m 2 t - u)Cf, 

92 = 2 Bq + 2B 7 q + (rrif — u)(2Cq + 3Cf) + 4(1 — e)Co 0 , 

g 3 = 2{e - 1 + Cl{ml - u) - (1 + e)B^ + 2 Bq + [-5Cf + 4Cf + 2(1 - e)^] 

+w[Cf — 4 C 2 — 2(1 — e)^^]}, 

g 4 = 4m t [Cl - 2C* - (1 - e)^]. (67) 
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For diagram (j) in Fig. 2, we can write the form factor as 


fi 


(C F - C a /2) 


87T 2 (m 2 — u ) 


(68) 


9i = m t (m 2 t - u) [Cl + e(Cl + C 2 7 )], 

g 2 = 2Clm 2 — (1 + e)Bl + + B 7 Q + (1 + e){m 2 — u)C'[ + 2(—1 + e)Cg 0 , 

(yf 3 = —1 — + 2 Sq + 2 rtCj 2 — 2m 2 (2Cl + C 22 ) 

—e[—1 + i?o — 2 (m 2 — w)(Cj + Cj 2 )], 

(74 = — 2m,t[Cle + (—1 + 2e)Cl + (—1 + c)Cl 2 ]. (69) 


For diagram (k) in Fig. 2, we can write the form factor as 


q 3 

fi = CF-^g, 
\b r K z s 


(70) 


m = f = (i - 

For diagram (1) in Fig. 2, we can write the form factor as 

fi = C F . ^ 2 - yZ~ 9 i 

lb7r z (rnf — u) z u 


(71) 


(72) 


fh = —m t (m 2 - u) {(e - 1 )m 2 t (Bl - Bf) + u(e - 3)B 7 } , 
92 = y = (~1 + e)mt(Bl - B 7 0 ) - (-1 + e)u 2 B 7 

+ m t M {(1 — e )[-^o — 2(1 + -BoD] + 2(—3 + e)Bg} . 

For diagram (m) in Fig. 2, we can write the form factor as 


(73) 


(74) 


9i — —m t [CQ + Cf + C 2 — 2(s — u + m 2 )(D q + Fh + _D 2 + -D 3 ) + 4_D 0 o], 
g 2 = -2(C' 0 4 + C 5 0 ) - 2D 3 0 (m 2 - u) + C{ + Cf + 2(t + 2m) (D 2 + D 3 ) 

~ m h±(Di + 2D 2 + D 3 ) — (Di + 4-D 2 + 3 -D 3 ) + 4e.Doo, 

S's = ~2 |2 [Cq — (1 + e)Co + Cg] — 2e(Cf + Cf) + Cf + Cf + 2.Df [— m 2 H ± + (1 — e)mf 
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+ (1 + e)w] + 2[(3 — 2e)m 2 — t + 2ne](Z)i + D 2 + -D 3 ) + 2 [(1 — e)m 2 + we] (.Du 

+ 2 D 12 + 2 D 13 + D 22 + 2 D 2 ,3 + -D 33 )} , 

9 a = 4m^ {(D > 2 + D 3 + -Du + 2 .Di 2 + -D 13 + -D 22 + -D 23 ) — e [-^o + -^11 + D 22 + -D 33 
+2(D>i + D 2 + -D 3 + -D 12 + -D 13 + -D 23 )]} , 

g 5 = 4 m t |Cq + 2eCi — C\ + 2(w — m^)D 1 + [m#± — (7 — 2 e)m 2 + 2 (f + 2 m) — 2we]-D 2 
+ [— Tn\j± ~ m 2 + 2 s]Z ?3 + [ 2 (e — l)m 2 — 2 we](.Di 2 + -D 22 + -D 23 ) j, 

9g = —4-m t [D >0 + D >1 + (2 - e)D 2 + D 3 + D 12 + -D 22 ], (75) 

where the variable of the D-function is the same with Dq. 

For diagram (n) in Fig. 2, we can write the form factor as 

fi = {C F - C A /2)-^ gi (76) 

9i = -mt {^[s + e(m 2 - w)] + s(D 2 + D 3 ) 

—e[C *2 — m^Di + D 3 ) + u(D 2 + D 3 )] — 2 .D 00 }, 
g 2 = — Cq — Oq + Cq (2 + e) + D^\m? H ± + em^ — (2 + e)w] + Di\jrn^{\ + e) 

+s — w(l + e)] + {D 2 + D 3 )[em 2 + s — w(l + e)] — 2eD 3 Q, 

93 — 2 j—C q — Cq — (Cj + Cj) — e(C{ + C®) + 2m 2 _D 2 — t(D 2 + D 3 ) + mf(2D 3 + 77 2 2 
+ 2 D 23 + -D 33 ) — e C 2 + (rrif — w) [D >2 + D 3 + D 22 + 2 D 23 + -D 33 ] } , 

9a = 2m* {— D 2 — D 3 — D\ 2 — D \ 3 — D 22 — 2D 23 — D 33 
+e[-Do + 2D ) 2 + 2D 3 + D 22 + 2D> 2 3 + -D 33 ]|, 

9b — —2 |(7g(—1 + e) + (4m 2 —t — 2u)D 2 — sD 3 + m%(D 22 + D 23 ) 

+e[ 2 C^ — + C® — (m 2 — u)(D 3 + D 2 2 + -D 2 . 3 )]} > 

96 = 2 m t [D 0 — (—3 + e)_D 2 + -D 12 — (—1 + e)(D 22 + -D 23 )], (77) 

where the variable of the D-function is the same with Dq. 

For diagram (o) in Fig. 2, we can write the form factor as 

/. = (C F - cy2)|7 9i (78) 

fid = —— cCq + -D^(l + e)s — e(C® + C®) + s(l + e)(Dh + D 2 + D 3 ) — 2D 0 o], 
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92 — ~C-q ~~ Co + Cq(2 + e) + Dq(— 2mf + t) — m 2 (D\ + D 2 ) 

+[s(l + e) — u\D 3 — 2eDoo> 

5-3 = 2 |Cq( 1 + e) - Cge + Dl[2m 2 H ± + m 2 (4 - e) - s(2 — e) - 3t - u{2 - e)] 

+e(Cj* + Cf — C ) 1 + C| ~ C 2 ) + (-Di + -D 2 + -D 3 ) [m^± (—3 + e) — ? 7 i 2 e + 3s + t (2 — e) 
+u (3 + e)] +(-Dn + 2 i 4 i 2 + 2Z>i 3 + D 22 + 2 .D 23 + D 33 ) [m 2 (l — e) + ew]|, 

54 = 2m t {— D 3 — Dis — D-23 — -D 33 + e(-D 0 + 2Z?i + 2D 2 + 2.D 3 + Du 
+2 /} 12 + 2 D 13 + D 2 2 + 2 D 23 + -D 33 )} , 

95 = —2 |Cq — -Dg(? 77 ^± — 2m 2 ) — C* + e(Cj ! + C[?) + -Di[—7(1 + e) — 2w] + 

(D 2 + -D 3 ) [em 2 ff ± — s(l + e)] + D :i [—t(l + e) — m] + (Du + _D 12 + 2.D 13 + D 23 + D 33 )ue 
+m 2 [5-Di + -D 2 + 4i4 3 — (—1 + e)(Z?n + D \ 2 + 2D V3 + H 23 + -D 33 )]} , 

9 o = 2rn t |-Dq — (—2 + e)(Z?i + -D 3 ) + _D 13 + _D 33 — e(-D n + D u + 2.D 13 + H 23 + -D 33 )} » (79) 

where the variable of the D-function is the same with Dq. 

By decomposition, the loop integrals in above form-factors can be calculated by the 
limit number of scalar integrals Bq , Co and Do. The scalar integrals Co and Do are UV 
finite, however some of them contain soft and collinear divergences. Because the finite 
scalar integrals could be calculated by numerical method [19], only the divergent ones are 
presented explicitly in this paper. It should be noted that only real part of the integrals, 
which is relevant to our results, is given. 

D 0 and C 0 scalar integrates could be generally written as 

Dq = (4vr/i 2 ) e r(l + e)(— H—-—f do), 

e e 

Co = — ( 47 T/ 7 2 ) € r(l + e)(— H-f Co). 

e- e 

The coefficients d 2 , d\ and do of D q are 

^ 2 2 (m t 2 — t) ( m t 2 — u) ’ 

, = _I_ m t (m 2 — m 2 H ±) 

1 (m t 2 — t) ( m t 2 — u) ( m t 2 — t)(m t 2 — u) ’ 


(80) 

(81) 


do — 


2 ( m t 2 — t) (m t 2 — u) 


7T 


0 (m H ± -m t ) - 


2 log 2 m t + 3 log 2 ( m 2 — t) 
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+3 log (rnf — u) — 2 log {mj — t) log(t) — 2 log (m 2 — u) log(w) — 3 log^ \m 2 H ± — m 


mf 


m 2 H ±s 




-m,( t, I-». 

\mt\u — rn H ±) 


m t 


m t 2 — rriH± 2 


+ 2/i 2 


u — 


m H± \ 


m t 2 — t) 


+ 21 % 


m t 


+2U J t -P>±)-2 l , 2 ( 

V — u / V 


(t - m^± 2 )(-u - m^± 2 )\ + / me 


{"'I ~ t)(nif ~ u) 


nit — u 


‘ \m 2 (m 2 H ± — u)(m 2 H± — t) 


ts 


—2 li<2 


m 2 s 


■T 2/i 2 


(m 2 H ± -t)(t- mt) 
su \ 1 


, (m 2 — t ) (m 2 — u)J V (m^± — w) (m — m 2 ). 
The coefficients d 2 , di and do of T>q are 

3 


do — 


2 s (t — m t 2 ) ’ 


m t 2 — t 
_ 2/j f <(“ - ”?) 

’UKi-i) 


(83) 


(84) 


di = 


log 


m t (jn 2 H ± - m 2 ) 


s (t — m t 2 ) s(m t 2 — t) 


2 ’ 


(85) 


do = ——^2 {it 2 [1 - 0(m H ± - mt)} + 21og 2 (m,) + log 2 |m^± - m 2 | - log 2 (m 2 - t) 


-2 log(m|* - <) Mm? - *) + 4 log log-^ - ^(-^-5 

+ 2US'A-1) - 2 ZA ) + g<,(-gk-) +2i!2 ((%^) 

V rrif / \mf(m z H± -t)J \m z H ±-tJ \ mf— t J 

The coefficients c 2 , Ci and c 0 of Cq are 


( 86 ) 


02 =-, 

s 


(87) 


Cl = 


log(s) 


( 88 ) 


Co = 7 T- [4vr 2 -31og 2 (s) 
os L 


(89) 


The coefficients Ci and cq of Cg are 


Cl -.2 


m H± ~ t 


log 


/ 9 O' 

_ 

V - t 


(90) 
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= 2{mjj ± - t) {'° g2(m? “ t] - l0g2(m ‘ “ ' 4i) + 2lk LC H -^ ) ~ 2 ' !2 (r^W)} ’ (91) 

The coefficients C 2 ,ci and cq of Cfj are 


1 

C2 2 {mi-ty 

1 , rn t 

Cl “ (m? - t) ° g m? - t' 

c ° = 2imhr){- log(m<)+log2(m? - () “ 2,!2 (r=W)} ’ 

Cq, Cq and can be obtained from Cq, Cq and Dq by replacing t with u. 


(92) 


(93) 


(94) 
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FIG. 1. Feynman diagrams at LO for bg —»■ tH 
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FIG. 
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. Feynman diagrams of the virtual correction for the process bg 


tH~. 
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FIG. 3 



. Feynman diagrams of gluon-radiation process of bg —► tH g. 
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FIG. 5. (a): the total cross sections as a function of fi/no with /io = m H ± + where tan/3 = 2 
and m H ± = 200 GeV. Curves (l)-(4) represent the cross section at NLO in OS scheme (1), NLO 
in MS scheme (2), LO in OS scheme (3), LO in MS scheme (4). (b): 6 (defined in text) as a 
function of ///Vo- 







Charged Higgs Mass (GeV) 


FIG. 6. Various contributions to K-factor versus m H ± with // = n 0 . Curves (l)-(6) represent 
contributions from all (1), improved Born (2), virtual plus gluon-radiation (3), qq (4), bq(q ) (5) 
and initial-gluon (6). 







